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Abstract. The remarkable efficiency of the UVES spectrograph at the VLT has made 
it possible to push high-resolution, high-S/N ground observations of the Ly-a forest 
down to z ~ 1.5, gaining new insight into the physical conditions of the intergalactic 
medium and its evolution over more than 90% of the cosmic time. The universal ex¬ 
pansion, the UV ionizing background and the gravitational condensation of structures 
are the driving factors shaping the number density and the column density distribu¬ 
tion of the absorbers. A (limited) contribution of UV photons produced by galaxies 
is found to be important to reproduce the observed evolutionary pattern at very high 
and low redshift. The Lyman forest contains most of the baryons, at least at 2 > 1.5, 
and acts as a reservoir for galaxy formation. The typical Doppler parameter at a fixed 
column density is measured to slightly increase with decreasing redshift, but the in¬ 
ferred temperature at the mean density is increasing with redshift. The signatures of 
Hell reionization and feedback from the formation of galactic structures have possibly 
been detected in the Lyman forest. 


1 The Observations 

A sample of 8 QSOs with 1.7 < Zem < 3.7 has been observed with 

VLT/UVES at a typical resolution 45000 and S/N ~ 40 — 50. Thanks to the 
two-arm design of the spectrograph Q, a high efficiency has been achieved in 
the whole optical range, from the atmospheric cutoff to 1/im, which translates 
immediately in the possibility of obtaining new results on the Lyman forest, 
especially at z ^ 2.5. The data have been reduced with the UVES pipeline Q - 
an non-negligible factor in maximizing the scientific output per unit time - and 
analyzed with the VPFIT package . 

2 The number density and column density evolution of 
Ly-CK lines 

The swift increase of the number of absorptions (and the average opacity) with 
increasing redshift is the most impressive property of the Ly-a forest. Fig. ^ 
shows the number density evolution of the Ly-a lines 00 in the column 
density interval 0 Nni = cm The maximum-likelihood fit to the 

^ This range in Nui has been chosen to allow a comparison with the HST Key- 
Programme sample at z < 1.5 for which a threshold in equivalent width of 
0.24 A was adopted. 
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Fig. 1. Number density evolution of the Lya forest with Nhi = cm“^. 

Dotted lines refer to the evolution compatible with an ionising UV background due 
only to QSOs. Solid lines show the expected evolution when both QSOs and galaxies 
contribute to the background, for models with /esc=0.05 (upper line), 0.1 and 0.4 (lower 
line). Data points come from several observations in the literature, as given by | [li| ]. 
The modelled evolution has been normalized to the observed evolution in the redshift 
range 2 < 2 < 3. 0 


data at z > 1.5 with the customary power-law parameterization provides N(z) = 
Noil + = (6.5 ± 3.8) (1 -b The UVES observations imply that the 

turn-off in the evolution does occur at z ~ 1, not at z ~ 2 as previously suggested. 

While the opacity is varying so fast, the column density distribution stays 
almost unchanged. The differential density distribution function measured by 
UVES |^,|^ , that is the number of lines per unit redshift path and per unit Nhi 
as a function of Nhi, basically follows a power-law /(Nhi) oc extending 

over 10 orders of magnitude with little, but significant deviations: the slope 
/3 of the power-law in the range 14 ^ log Nhi ~ 16 goes from about —1.5 at 
< z >= 3.75 to —1.7 at z < 2.4. Recent HST STIS data ^ confirm that this 
trend continues at lower redshift, measuring a jS of —2.0 at z < 0.3. 

3 The evolution of the Ly-a forest and the ionizing 
background 

The evolution of the N{z) is governed by two main factors: the Hubble expansion 
and the metagalactic UV background (UVB). At high z both the expansion, 
which decreases the density and tends to increase the ionization, and the UVB, 
which is increasing or non-decreasing with decreasing redshift, work in the same 
direction and cause a steep evolution of the number of lines. At low z, the 
UVB starts to decrease with decreasing redshift, due to the reduced number 
and intensity of the ionizing sources, counteracting the Hubble expansion. As a 
result the evolution of the number of lines slows down. 
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Up to date numerical simulations [ p^ have been remarkably successful in 
qualitatively reproducing the observed evolution, however they predict the break 
in the dN/dz power-law at a redshift 0 ~ 1.8 that appears too high in the 
light of the new UVES results. This suggests that the UVB implemented in the 
simulations may not be the correct one: it was thought that at low redshift QSOs 
are the main source of ionizing photons, and, since their space density drops 
below z ~ 2, so does the UVB. However, galaxies can produce a conspicuous 
ionizing flux too, possibly more significant than it was thought |]^. The galaxy 
contribution can keep the UVB relatively high until at z ~ 1 the global star 
formation rate in the Universe quickly decreases, determining the qualitative 
change in the number density of lines. 

Under relatively general assumptions, it is possible to relate the observed 
number of lines above a given threshold in column density or equivalent width 
to the expansion, the UVB, the distribution in column density of the absorbers 
and the cosmology |^: 

= C[{l + zfr-l{z)\^-^ H-\z), (1) 

>Nhi ,lim 

where /hi is the photoionization rate and P the power-law index of the Nhi 
distribution. 

To estimate /hi we have investigated the contribution of galaxies to the 
uvb|), exploring three values for the fraction of ionizing photons that can es¬ 
cape the galaxy interstellar medium, /esc = 0.05,0.1 and 0.4 (the latter value 
corresponds to the Lyman-continuum flux detected by [ p^ in the composite 
spectrum of 29 Lyman-break galaxies). Measurements of the UVB based on the 
proximity effect at high-z and on the Ha emission in high-latitude galactic clouds 
at low-z provide an upper limit on fesc ~ 0.1, consistent with recent results on 
individual galaxies both at low-z and at z 3 [|j. Introducing a con¬ 

tribution of galaxies to the UVB, the break in the Ly-a dN/dz can be better 
reproduced than with a pure QSO contribution . The agreement improves con¬ 
siderably also at z ^ 3. Besides, models with Qa = 0.7, Qm = 0.3 describe the 
flat evolution of the absorbers much better than Qm = 1. A consistency check 
is provided by the evolution of the lower column density lines. For XogNui ~ 14 
the Nhi distribution follows a flatter slope /3, and according to Eq. this trans¬ 
lates directly into a slower evolutionary rate, which is consistent with the UVES 
observations]!^: dN/dz(i 3 i<jvh/<i4) oc (1 -I- Another diagnostic can 

be derived from the spectral shape of the UVB and its influence on the intensity 
ratios of metal lines 00 - 

4 Mapping the column density distribution into the mass 
distribution of the gas 

It is instructive to transform the observed column density distribution in the 
mass distribution of the photoionized gas (Fig. and interpret it, following 
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Fig. 2. The differential mass density distribution of the Ly-a forest as a function of 
Nhi- The arrow indicates the direction towards which the points move if the temper¬ 
ature or the ionization rate increase 


Schaye|^, as a function of the matter density contrast: 1) the flattening at 
log Nhi ~ 13.5 is partly due to line crowding and partly to the turnover of the 
density distribution below the mean density; 2) the steepening at log Nhi ~ 14, 
with a deficiency of lines that becomes more and more evident at lower z, reflects 
the fall-off in the density distribution due to the onset of rapid, non-linear col¬ 
lapse 3) the flattening at Nhi ~ 10^® cm“^ can be attributed to the flattening 
of the density distribution at density contrast ^ 10^ due to the virialization of 
collapsed matter. The differential mass density distribution has a sort of uni¬ 
versal form when plotted as a function of the density contrast. A given density 
contrast, however, corresponds to lower and lower column densities with de¬ 
creasing redshift, and this causes a translation of the mass density distribution 
(Fig. 1) towards the left with decreasing redshift, which explains the steepening 
of the slope [3 reported in Sect.2. Hydrodynamical simulations successfully re¬ 
produce this behaviour, indicating that the derived matter distribution is indeed 
consistent with what would be expected from gravitational instability. 
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5 The cosmic baryon density 

The amount of baryons required in a given cosmological scenario to produce 
the observed opacity of the Lyman forest can be computed Q under general 
assumptions. A lower-bound to the cosmic baryon density can be derived from 
the mean Ly-a flux decrement, D,|p^ and/or from the distribution of the Ly-a 
optical depths. The limits derived from the effective optical depths measured in 
the UVES spectra at 1.5 < z < 4 are reported in Tab. |^. They are consistent 
with the BBN value for a low D/H primordial abundance. Most of the baryons 
reside in the Lyman forest at 1.5 < z < 4 with little change in the contribution 
to 12 as a function of z. Conversely, given the observed opacity, a higher UVB 
requires a higher 12;,. As pointed out by |P, an escape fraction as large as 0.4, 
as measured by would result in an 12b ~ 0.06 in conflict either with the 

primordial D/H abundance or in general with the BBN or with the Ly-a opacity 
measurements. 


Table 1. Lower limits to 12b/i^'® derived from the UVES spectra at 1.5 < z < 4 (for a 
universe with 12m = 0.3,12yi = 0.7) 


UVB 

T = 2 ■ 10'‘ K 

T = 6 • 10^ K 

QSOs 

0.017 

0.011 

QSOs -f GALs 

(fe.c = 0.1) 

0.028 

0.018 


6 The temperature of the IGM 

If the Ly-a forest is in thermal equilibrium with the metagalactic UV back¬ 
ground, the line width of the absorption lines, described by the b parameter 
of the Voigt profile, is directly related to the gas temperature of the absorbing 
medium determined by the balance between adiabatic cooling and photo-heating: 
b = \/2kT jraion- Additional sources of broadening exist, such as the differential 
Hubble flow across the absorbers, peculiar motions, Jeans smoothing. However, 
there is a lower limit to the line widths, set by the temperature of the gas, that 
is in principle measurable. In practice the situation is slightly more complex 
because for a photoionized gas there is a temperature-density relation, the so- 
called equation of state: T = To (1 -I- where T is the gas temperature. 

To is the gas temperature at the mean gas density, Sb is the baryon over-density, 
(pb — 'Pb)l'Pb 7t is a constant which depends on the ionization history. The 
equation of state translates into a lower cutoff bciNni) in the Nni-b distribu¬ 
tion. 

The observed cut-off Doppler parameter at a hxed column density of logjy^^ = 
13.6, 6c(13.6) is measured to increase with decreasing redshift, while the slope 
of the cutoff does not change significantly. A typical value of be ~18 km s“^ at 
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1.5 < z < 4 corresponds to a reference temperature of 2 • 10‘^K. This does not 
mean that the temperature of the IGM increases with decreasing z: on the con¬ 
trary, taking into account the equation of state and the fact that a given column 
density corresponds to higher and higher over-densities with decreasing redshift, 
it turns out that the temperature at the mean density is actually decreasing with 
decreasing z. Furthermore, evidence has been found p^ | for a general increase 
of the temperature around the redshift z = 3.3 ± 0.15, attributed to the reion¬ 
ization of Hell. Temperature and ionization fluctuations are also expected due 
to feedback processes from the formation of galactic structures |]^ and might 
have been observed in correspondence of some voids in the Ly-a forest []ll|. 
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